Bubble domains in disc-shaped ferromagnetic particles 
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We study the fundamental magnetic states of disc-shaped ferromagnetic particles with a uniaxial 
anisotropy along the symmetry axis. Besides the monodomain, a bidomain state is also identified 
and studied both numerically and theoretically. This bidomain state consists of two coaxial oppo- 
sitely magnetized cylindrically symmetric domains and remains stable even at zero bias field, unlike 
magnetic bubbles in ferromagnetic films. For a given disc thickness we find the critical radius above 
which the magnetization configuration falls into the bidomain bubble state. The critical radius 
depends strongly on the film thickness especially for ultrathin films. The effect of an external field 
is also studied and the bidomain state is found to remain stable over a range of field strengths. 
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Current developments in experimental techniques al- 
low the engineering of magnetic particles to very small 
length scales, typically tens or a few hundred of nanome- 
tres pj. In such systems the characteristic length scales 
of the magnetic domains may be comparable to the sys- 
tem size. The particles should thus support only a small 
number of magnetic domains and these should be usu- 
ally biased by the symmetry of the particle to give a 
small number of distinct magnetic configurations. Exam- 
ples are vortex states observed in disc particles 0, 0, Q 
as well as vortex and onion states observed in magnetic 
rings 0,13- This situation is in marked contrast to a bulk 
ferromagnet where we have a large number of magnetic 
domains with significantly different sizes. 

Recently, much research effort has focused on the inves- 
tigation of narrow ferromagnetic rings where it is found 
that single, highly symmetric domain structures mediate 
the switching between the high moment (onion) state and 
low moment (vortex) states Motivated by the work 
on rings, in the present case we search for the existence 
of high symmetry stable states in circular disc structures 
with perpendicular anisotropy, in view of the recent de- 
velopment of materials with a large uniaxial anisotropy 
0, 13 and experiments in permalloy discs [t| . In partic- 
ular, we identify bidomain states which exist even in the 
absence of a bias field but are otherwise the analogues of 
the magnetic bubbles observed in ferromagnetic films. In 
addition, the bidomain states remain stable for a range 
of applied fields and can be manipulated in a controlled 
manner. Furthermore, we estimate the size of the mag- 
netic domains supported in the particle. The importance 
of the present findings is that, again, very simple, high 
symmetry domain structures are found to be stable and 
to mediate the magnetic switching process, in contrast 
to the complex behaviour which usually occurs in small 
elements. 

Static as well as dynamical properties of the magneti- 
zation m are governed by the Landau-Lifshitz equation. 
The constant length of the magnetization is normalized 



to unity: m 2 — 1. An important length scale of the 
system is the exchange length 
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where A is the exchange constant and Mq is the satu- 
ration magnetization. In the following, we shall use £ ex 
as the unit of length. Another important quantity is the 
dimensionless quality factor 
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where K is the anisotropy constant. The significance of 
the quality factor can be seen in two important quanti- 
ties. First, the domain wall width is £ cx / y/n. Second, k 
controls the relative strength of the magnetostatic field 
which has a demagnetizing effect, with respect to the 
anisotropy field which favours alignment along a direc- 
tion perpendicular to the film. We shall suppose here 
that k > 1 which means that the anisotropy is, in gen- 
eral, stronger than the demagnetising field. 

The Landau-Lifshitz equation is the basis for all our 
calculations and we are interested only in its static so- 
lutions. We find such magnetic configurations by a re- 
laxation algorithm the details of which were explained 
in ^(J |n|. We introduce a Gilbert damping term in 
the equation and feed the algorithm with an initial guess 
state. This eventually converges to a static solution at a 
local minimum of the energy functional. 

The most demanding part of the method is by far 
the calculation of the magnetostatic field which requires 
the solution of a boundary value Poisson problem. This 
makes the numerical simulation of most realistic prob- 
lems practically impossible to achieve in three dimen- 
sions. However, a huge reduction of the numerical cal- 
culations is obtained if we confine our interest to axi- 
ally symmetric configurations, and in the rest of this 
paper we shall be concerned only with such magnetic 
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states. We expect that this is not a serious constraint 
at least for the most basic magnetic states of small par- 
ticles which will be our main focus. Indeed, it is rea- 
sonable to assume that the lowest lying states of the 
system will have the symmetry of the geometry of the 
particle. We call z the axis of symmetry of the disc 
which is also the direction of the easy axis. We go to 
cylindrical coordinates and suppose that the radial (m p ) , 
azimuthal (m^) and longitudinal (m z ) components of 
the magnetization vector are functions of p and z only: 
m p = m p (p, z), m = m^(p, z), m z = m z (p, z). 

As a first step we calculate the simplest possible state. 
This is expected to be a single-domain state in which all 
spins are driven by anisotropy and lie roughly along the 
symmetry z axis. We use the uniform m = (0, 0, 1) state 
as an initial guess in the relaxation algorithm. This then 
quickly converges to a quasi-uniform static state, at least 
for strong anisotropy K > 1. The magnetization vector 
deviates from the z axis only around the edges of the 
disc. The quasi-uniform state is a monodomain state and 
is thus expected to be the ground state for sufficiently 
small particles. The transformation m — > — m gives a 
second monodomain state. 

As the size of the particle becomes larger it is antici- 
pated that the magnetic configuration will break up into 
domains. We conjecture a bidomain state that is axi- 
ally symmetric. This consists of an inner cylindrical do- 
main of "down" magnetization surrounded by the outer 
domain of "up" magnetization. A domain wall has to 
separate the two domains. We shall call these axially 
symmetric bidomain states "bubbles" because they bear 
some essential similarities to the so-called magnetic bub- 
bles observed in abundance in ferromagnetic films. 

Ferromagnetic films with a strong perpendicular 
anisotropy were studied experimentally and theoretically 
around the 70s. These early studies were largely driven 
by technological interest in magnetic bubbles whose stat- 
ics and dynamics were analysed in detail |12| . A magnetic 
bubble is a circular domain of opposite magnetization in 
an otherwise uniformly magnetized film with magnetiza- 
tion perpendicular to the film. The presence of an ex- 
ternal bias field is essential for the stabilisation of these 
structures 0, Q . If the bias field is lifted then the mag- 
netostatic field destroys the bubble which expands and 
eventually transforms into stripe domains. In contrast, 
the bubble states calculated here for sufficiently small fer- 
romagnetic particles remain stable even in the absence of 
a bias field. 

We thus return to our conjecture and proceed to test 
it numerically. As a standard example we choose the 
following values for the quality factor, the film thickness, 
and the radius: 

K = 2, d = 8£ eX) R = 2A£ ex , (3) 

and this choice will be explained later in the text. In 
the case of the FePt films of Ref. 0] where M = 
1150 emu/cm 3 , A = 10~ 6 erg/cm, the exchange length is 
t ex = 3.5 nm and thus the values of Eq. J3J are translated 
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FIG. 1: The bubble state illustrated at the middle (z — 0) 
plane of the disc and at the top (z = d/2). The magnetization 
is such that m z = — 1 at the center and m z m 1 in the outer 
domain. The arrows show the projection of the magnetization 
on the (x, y) plane which has a significant value at the domain 
wall. This is a Bloch-like wall in the middle (z = 0) plane 
and it turns almost Neel-like at the top and bottom surfaces 
(z = id/2). The magnetization satisfies the parity relations 
m p (p,z) = -m p (p,-z), m<p(p,z) = m#(p, -z), m z (p, z) = 
m z (p, —z). 



to d = 28 nm, R — 84 nm. The value n — 2 for the qual- 
ity factor corresponds to K = 1.6 x 10 7 erg/cm 3 for the 
anisotropy constant, as is typical in FePt films. We feed 
our numerical algorithm with an initial condition which 
has the gross features of the bubble state described above 
with a domain wall of the Bloch type smoothly connect- 
ing the domains. The algorithm converges to a static 
bubble state which has a complicated domain wall struc- 
ture shown in Fig. ^ Also, the magnetization deviates 
to some extent from the z direction at the side surface of 
the particle. The profile of this structure is sufficiently 
interesting and deserves attention. The magnetostatic 
energy is the driving force here and it clearly favours a 
bidomain state with opposite magnetization where the 
total magnetization would roughly vanish. On the other 
hand, the anisotropy and exchange energies are signifi- 
cant at the domain wall and they put a tension on it to 
shrink. In the final result the bubble has an inner do- 
main with volume smaller than the outer domain. Thus 
the total magnetization points along the symmetry axis 
and it is nonzero. As mentioned already, contrary to the 
situation in films, we suppose here that no external bias 
field is present. 

The domain wall resembles those discussed in the lit- 
erature in related calculations 0, 0|. It is Bloch in 
the central plane (set at z = here) and it progressively 
becomes Neel towards the surfaces. The Neel wall is sig- 
nificantly wider than the Bloch wall. The radius of the 
bubble is larger at the centre than near the surfaces, but 
this is a small effect. It is easy to understand that for 
this type of wall the magnetostatic energy and the total 
energy density are larger near the surfaces than at the 
disc centre. In short, the surfaces disfavour the domain 
wall and subsequently also the bubble state. The final 
and important result is that, for the parameters @, the 
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FIG. 2: Energy per unit volume (in units of 2ttMq) of the 
monodomain (dashed lines) and of the bubble state (solid 
lines) as a function of the disc radius R for three values of 
the disc thickness d (R and d in exchange length units). The 
bubble exists only for R greater than a critical radius R\ and 
it has a lower energy for R > R c where R c is yet another 
critical radius which corresponds to the intersection of the 
two lines for each value of d. Both R\ and R c depend on d. 



bubble state has a lower energy than the monodomain 
state. 

We now proceed to a systematic numerical study of the 
bubble state. We first fix the disc thickness at d — 8.0 
and vary the radius R. We find a bubble state when 
the radius is larger than some critical radius R\ m 14. 
For smaller radii R < R\ our algorithm always converges 
to the quasi-uniform state irrespective of the initial con- 
dition. Our results thus indicate that the inner bubble 
domain cannot be sustained if it is too small. For a disc 
radius slightly larger than Ri the bubble radius is small 
and the total magnetization of the structure is large. As 
the radius of the disc increases the inner bubble domain 
expands and the absolute value of the total magnetiza- 
tion decreases. 

The energy per unit volume of the bubble state as a 
function of the disc radius is given in Fig. [21 for three 
values of the thickness d, along with the corresponding 
energy for the monodomain state. The latter exists as 
a local minimum of the energy for any radius up to the 
largest that we checked. The energy per unit volume 
of the bubble is greater than that of the monodomain 
state at the lowest radius R\ where the bubble first ap- 
pears. It decreases for larger radii and becomes lower 
than that of the monodomain state above a critical ra- 
dius R c . It eventually becomes an increasing function of 
the radius but apparently remains lower than the energy 
of the monodomain state for all R > R c . Wc consider R c 
as marking the size for the break up of the magnetization 
configuration into domains. We could find the bubble as 
a local minimum of the energy for all radii R > R\ that 
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FIG. 3: The critical radius R c as a function of the thickness 
of the disc d. For R > R c the bubble has a lower energy than 
the monodomain state. 

we have checked. However, it is expected that a multi- 
domain state will eventually set in for a sufficiently large 
radius, with energy lower than the energy of both the 
monodomain and the bubble. 

In order to study the dependence of R c on the thickness 
d we have repeated our calculation for a few values of d. 
From Fig. [5] one can extract the critical radii for three 
values of d = 6, 8, 15. In Fig. |21 we give the R c as a 
function of d inferred from six values of d. For small d the 
critical radius significantly exceeds the particle thickness 
because surface effects become important and disfavour 
the formation of a domain wall. On the other hand, the 
critical radius R c levels off for higher values of d. The 
thickness for which R c attains a minimum appears to be 
close to d = 8 for which R c as 18. This is actually the 
reason for choosing @ as our standard parameters along 
with the fact that for d = 8 the energy of the bubble has 
a minimum at R as 24 as is seen in Fig. 

We have also repeated our calculation for the particle 
sizes employed in the experiment of Ref. |9| and have 
confirmed the existence of a bubble state. However, a 
detailed quantitative comparison cannot be made before 
one determines the strength of the deposition induced 
anisotropy in the permalloy used in the experiment. 

Once we have established the existence of a bubble we 
would like to know how it behaves under an externally 
applied field. Apart from the apparent practical impli- 
cations, this is an interesting question also because the 
field will affect the intricate balance of energies that is 
responsible for the stabilization of the bubble. The field 
is applied along the symmetry axis of the disc, i.e., it is 
of the form h cxt = (0, 0, h cxt )- 

We apply the field on a particle which is already in 
a bubble state. As a specific example, we choose our 
standard parameter values ©. Our results are given in 
Fig. For h cxt = the total magnetization per unit 
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FIG. 4: The total magnetization fi per unit volume as a 
function of the applied field, for a bubble state in a disc of 
thickness d = 8 and radius R = 24. 

volume, /i = 1/V J m z dV, is non-zero. If we choose the 
inner domain to point down then [i is positive (Point 
C in Fig. UJ. Applying a positive external field favours 
the outer domain, which expands at the expense of the 
inner domain. The system does reach a new equilibrium 
state which is again a bubble state but with a smaller 
radius. This corresponds to an increased value for \i. For 
a high enough magnetic field the bubble becomes too 
small and it cannot be sustained by the system. In our 
example fx jumps to unity for /i cxt > hs — 0.13 which 
signals that the bubble shrinks to zero radius resulting 
in a monodomain state with magnetization pointing up 
(Point B). On the other hand, if we start from point C but 
now reduce h cxt to negative values the inner domain is 



favoured and pushes the domain wall to a larger radius, 
which is reflected in a smaller value for /j,. The total 
magnetization fi crosses zero for h cxt = —0.08, it then 
becomes negative and eventually fi jumps to minus unity 
for h cxt < Ha — —0.26 which means that the system 
is in the monodomain state pointing down (Point A). 
Below Ka the domain wall is attracted by the side surface 
and is expelled from the disc. The dashed line in Fig. 0] 
corresponds to the equivalent physical situation obtained 
by the symmetry transformation m — ► —m. h cxt — > 

The bubble is stable in the range Iia < h cxt < hs 
in which reversible behaviour occurs. For h < Iia and 
h > Kb irreversible jumps in the magnetization occur 
which correspond to the domain wall being attracted to 
the edge (point A) or shrink to the center of the disc 
(point B). The size of the jumps in the magnetization re- 
flect the size of the domain wall being anihilated (larger 
when the inner domain expands) and constitutes yet an- 
other example of how the geometry of the element con- 
strains the shape of the domain wall and the details of 
the switching process 0,0]. On the other hand, if a 
particle with R > R c is saturated by a strong in-plane 
field, it will eventually relax into a bubble state after the 
field is removed 0. 

In conclusion, we have studied the fundamental 
states of disc-shaped magnetic particles with uniaxial 
anisotropy along the axis of the disc. A magnetic bubble 
state has been identified within our numerical calculation 
and has been studied in detail. The bubble is a particu- 
larly simple axially symmetric state and has apparently 
been observed experimentally Q. 

We are grateful to N.R. Cooper, A. Ntatsis, CA. Ross 
and T. Shima for discussions. This work was supported 
by EPSRC grant no. GR/R96026/01 (SK). 



[1] C.L. Dennis, R.P. Borges, L.D. Buda, U. Ebels, J.F. 

Gregg, M. Hehn, E. Jouguelet, K. Ounadjela, I. Petej, 

I.L. Prejbeanu and M.J. Thornton, J. Phys. C: Condens. 

Matter 14, R1175 (2002). 
[2] T. Shinjo, T. Okuno, R. Hassdorf, K. Shigeto, and T. 

Ono, Science 289, 930 (2000). 
[3] A. Wachowiak, J. Wiebe, M. Bode, O. Pietzsch, M. Mor- 

genstern, and R. Wiesendanger, Science 298, 577 (2002). 
[4] M. Klaui, C.A.F. Vaz, L. Lopez-Diaz and J.A.C. Bland, 

J. Phys.: Condens. Matter, 15, R985 (2003). 
[5] F.J. Castano et al, Phys. Rev. B 69, 144421 (2004). 
[6] J. Rothman, M. Klaui, L. Lopez-Diaz, C.A.F. Vaz, A. 

Bleloch, J.A.C. Bland, Z. Cui, R. Speaks, Phys. Rev. 

Lett. 86, 1098 (2001). 
[7] T. Shima, K. Takanashi, Y.K. Takahashi, K. Hono, G.Q. 

Li, S. Ischio, J. Magn. Magn. Mat. 266, 171 (2003); 

G.Q. Li et al, J. Appl. Phys. 94, 5672 (2003). 
[8] J.A. Christodoulides, Y. Huang, Y. Zhang, G.C. Haji- 

panayis, I. Panagiotopoulos, D. Niarchos, J. Appl. Phys. 

87, 6938 (2000). 



[9] P. Eames and E. Dan Dahlberg, J. Appl. Phys. 91, 7986 
(2002). 

[10] S. Komineas, N. Papanicolaou, Physica D 99, 81 (1996). 

[11] S. Komineas, Ph.D. Thesis, University of Crete. 

[12] A. P. Malozemoff and J.C. Slonczewski, "Magnetic do- 
main walls in bubble materials" (Academic Press, New 
York, 1979). 

[13] A. A. Thiele, Bell Syst. Tech. J. 48, 3287 (1969). 

[14] A. A. Thiele, Phys. Rev. Lett. 30, 230 (1973); J. Appl. 

Phys. 45, 377 (1974). 
[15] T.G.W. Blake and E. Delia Torre, J. Appl. Phys. 50, 

2192 (1979); E. Delia Torre, C. Hagediis, and G. Kadar, 

AIP Conf. Proc. 29, 89 (1975). 
[16] C.A.F. Vaz, L. Lopez-Diaz, M. Klaui, J.A.C. Bland, 

T.L. Monchesky, J. Unguris, Z. Cui, Phys. Rev. B 67, 

140405 (R) (2003). 
[17] M. Klaui, C. A. F. Vaz, J. A. C. Bland, T. L. Monchesky, 

J. Unguris, E. Bauer, S. Cherifi, S. Heun, A. Locatelli, L. 

J. Heyderman, Z. Cui, Phys. Rev. B 68, 134426 (2003). 



